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TRAESLATIOML MOTION OF BODIES UNDER THE FREE SURFACE OF A 
HEAVY FLUID OF FINITE DEPTH* 

By M. D. Haskind 


In reference entitled "The Two-Dimensional Problem of the 
Vibration of Bodies under the Surface of a Heavy Fluid of Finite Depths" 
the prpblem was to determine the wave motion of a heavy fluid excited 
by the periodic vibrations of a body of arbitrary shape situated under 
the free surface of the fluid of finite depthj the method of N. E. Kochin 
(reference 2) was used. 

In the present paper^ the two-dimensional problem of the wave 
motion produced in a heavy fluid of finite depth by the horizontal 
rectilinear and uniform motion of a solid body of arbitrary shape 
immersed under the sirrface of the fluid is considered by the same 
method. 


1. Statement of the Problem 

The problem of the translatory motion of a solid body uinder the 
free surface of a heavy inconpressible fluid of finite depth will be 
considered. The case in which the motion of the body occurs with con- 
stant horizontal velocity c will be studied. The motion of the fluid 
will be defined with reference to a moving system of coordinates Oxy 
fixed to the body^ the x-axis coinciding with the undisturbed level of 
the fluid and directed along the direction of motion of the body^ and 
the y-axls directed vertically upward. 

It will be assumed that the motion of the fluid is potential and 
steady relative to the body. From the integral of Lagrange for the 
pressure within the fluids 

p - PO = pc - p v2 - pgy (l.l) 
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where Pq is the atmospheric pressure^ p the density of the fluids 

g the acceleration of gravity: <P(x,y) the potential of the absolute 
motion of the fluid, and v = [ grad «p| the magnitude of the absolute 
velocity of the fluid. 

The function cp(x,y) is determined from the bomdary conditions; 
the flow condition on the wetted contour of the body, 

•^ = c cos(n,x) on C (l«2) 

where n is the outer normal to the contour C; 
on the free boundary p = Pq, and hence 

^ ll - 1 - sy = ^ 

on the bottom of the channel for y = -hQ, the following condition 
applies 

= 0 (1.4) 

According to the theory of waves of small anplitude, condi- 
tion ( 1 . 3 ) may be linearized. For this purpose the boundary condi- 
tion ( 1 . 3 ) is referred to the x-axis and the term v^/2 neglected. In 
place of condition (l.3), 

|| - f y(x) = 0 (1.5) 

It is easily seen that on the free surface the following relation 

holds 

cy(x) = cp + const (l»6) 

where ilr is the stream function. In fact, when the stream function of 

the motion of the fluid relative to the body is denoted by ijfQ, there 
is obtained 




to = t - cy 


cy = t - to 


or 
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From this relation;, equation (l.6) follows, since the boundary of 
the fluid in the relative motion is represented by stream lines on which 
ijfQ is constant. For the free surface, it may be assijmed that \jfQ = 0. 

Hence, on the free surface. 


cy(x) = t 

and therefore boundary condition (1.5) assumes the form 


where 


S<P 


for y = 0 


(1.7) 



( 1 . 8 ) 


From condition (l.S) it is seen that the equation of the free sur- 
face will be 


y(x) 


g 



( 1 . 9 ) 


2. Fundamental Fonnulas of the Problem 

The problem may be mathematically formulated as follows. It is 
required to determine the characteristic function w(z) = «p + ii|r 
(z = X + iy; 1 = a/:T) , satisfying the conditions: 

1. For 0>y>-hQ in the region occupied by the fluid, the deriv- 
ative dw/dz is finite and at infinity for x -► +» , the derivative 
dw/dz vanishes. 

2. On the contoirr C, the smooth flow condition applies 

^ = c cos(n,x) 

3. On the free surface for y = 0, the linearized condition holds 
with regard to the constancy of the pressure 

Re (dw/dz + ivw) = 0 
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4, On the hottom of the channel for y = -hQ^ the folio-wing condi- 
tion holds 

Im dw/dz = 0 

In the region occupied hy the fluid, the point z is taken and 
two contoxars and are drawn, of which €» contains both the 

point z a.rifl the contour C, while the eon-tains the contour C, 

but not the point z (fig. l) • By the formula of Cauchy for a single- 
-valued function dw/dz = v(z). 



( 2 . 1 ) 


where the bar over a letter indicates, as usual, the transition to the 
con^ilex conjugate value. The follo-wlng notation is introduced 




•v(0 

z - ^ 


( 2 . 2 ) 


It is evident that V 2 _(z) is a holomoipjhic function in the entire 


plane of the complex variable outside the contour Cq_, having at infin- 
ity the order z~l and capable of being continued analytically in the 
entire part of the complex variable plane which lies outside the con- 
tour C, while V 2 (z) is a holomorphic function within the contour C*, 

by the extension of which an analytical continuation of this fxmetion 
may be obtained over the entire strip 0>y> -hQ. 


The function ¥ 2 ( 2 ) inay be represented in another form. For this 
it is possible to find a f-unction (a(z), which in the strip 0>y>-hQ 

has a single pole of the first order ^ + iT) -vd.th residue A/2jti 

and which satisfies conditions 1, 3, and 4. 

In fact, for a vortex of strength T , located at the complex point 
^ = g + iT], an expression for the complex velocity was obtained by 
Tikhonov (reference 3) 
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P P 

2jti(z - “ 2jri(z - ^ + 2ihoj “ 

sh X(ti + ho) cos X(z - g + iho) 
+ X) e^(- Xho) v^Xho ' - X ' chX^ 


sh Xo(t) + ho) , ^ 

Tv — ^ r-JT-T-r— sin q(z ~ ^ + iho) 


^ vho - ch2 Xolio 
■where \q is the real and positive root of the eq-uation 

V sh Xho = Xch Xho 


(2.3) 


(2.4) 


For c^< in all cases -where the function to he integrated 

has a singulai-ity, the principal value in the sense of Cauchy is taken 
under the integral. 


For c^ > equation (2.4) has only imaginary roots and the 

fourth term of formula (2.3), which determines the presence of free 
■waves, is absent. 


For a source of strength Q located at the complex point 
^ ^ + If], the expression of the complex velocity may be obtained in 

the same manner as in the case of a vortex. Without the computations, 
the final res^llt is 


2Jt(z - i] 2n{z - i + 2iho) ^ 



+ X } exp(- Xho) 


ch X(il + ho) sin \{z - % + iho) 
V sh Xho - X ch Xho 


dX - 


ch Xo(il + ho) 

Qv cos Xo(z - C + iho) (2.5) 

Vho - 


By the use of expressions (2.3) and (2.5), to obtain the function 
cd(z) may be obtained -without difficulty. For this purpose, since 
A = r + iQ, the follpwing expression is obtained after sinple 
transformations : 
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( v + 


f \ _ A A 

~ 2jti(z - 23tl(z - ^ + 2itiQ) 

A sin X(z “ ^ + 2 Uiq) - A sin x(z - 
e3cp( XIiq) ^ __ 

^ A cos “ C + 211 iq) - A cos Xq(z - 


21 


vho “ 


cL\ + 


( 2 . 6 ) 


Here^ as in the preceding formulas, the fourth term, which deter- 
mines the presence of free waves, is present only if c^< 

¥hen A = v(^) d^ is substituted in the previous formula and 
integration is carried out over the contour Cp, 


y(z) = Vi(z) -2-^ 


S Xq(z - ^ + zihg) 

vlp - 


{/ 


Bln \(z ^ + 2ihQ) 

+ X) exp(- XIiq) — ^ , v:, dX 


V sh Xho - Xch XhQ 




( V +X) e3tp(- Xho) 


sin X(z - U 

V sh Xhfi - X ch Xhr 


dX - 


cos XqCz - £) . _ 
Jtv ;; Vd^ 


vho - oh^ Xohg 


:} 


( 2 . 7 ) 


^ f ,a> 


If both points z and ^ are situated in the strip 0>y> -hQ, 
the following equation holds 


- 5 + 211k,)] 


( 2 . 8 ) 


With this equation taken in account, it is foimd from equation (2.7) 
that the function ¥3(2) can be represented in the form 
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V2(z) 


_1 

2)tl 


^ 6351 j^iX(z - t + 2iliQ)J 


( V + \ ) exp(- Xho) 


sin X(z - ^ + 21ho) 


dX - jtv 


COB Xq(z - ^ + 2ihQ)~l _ 


2jti 


r{^ 



V ah XhQ - X ch XhQ 

(v + X) exp(- Xho) - 

^ V eh Xhn - X ch Xh, 


vho - ch^ Xoho 


?"]«• 


sin X(z - n ^ cos Xq(z - C) 

— i dX - iw ^ 

vho - ch^ XqUo 


:]« 


( 2 . 9 ) 


The conjugate coD5)lex functions are introduced for real X 


H(x) = / v(0 exp - iXCdC, H(x) = / v(0 exp iX^d^ 

Jci vJci (2.10) 

By an interchange in equation (2.9) of the order of integration, 
and by siiiiple transformations, there is readily obtained 



( V + X ) 635 (- Xho) 


(h(- X) 635 iX(z + 2iho)J - H(X) 035 1X(z + SlhQ^ - 


2(v sh Xhg - X ch Xho) 

H(X) 635 iXz + H(- X) 635(- iXz)^ 


dX - 


JtXV 


2(vh - oh^ 


( H(- Xq) 635 iXo(z + 2iho) + H(Xq) 635 - iXo(z + 2iho) - 


H(Xo) 635 iXQZ - H(- Xq) 635(- iXo^i) 


)] 


(2.11) 


It is of interest to find the character of the waves that remain 
behind the moving body. For this purpose the asymptotic expression of 
the complex velocity is first obtained for x-*-- a> in the case of a 
vortex and source. In reference 3, the asymptotic expression of the 
complex velocity in the case of a vortex is of the form 
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shXo(il + Iiq) , , / ^ 

« = - 2Fv P sinXo(z - ? + Iho) (2..12) 

^ “ vlio - ch^ XqIio 

In a similar manner, the asyniptotic expression of the complex 
velocity is obtained in the case of a source. Without the computations, 
the final result is 


ch Xo(n + bo) , ^ , 

(“q)„ ^ «= - 2Qv , p ' cosXqCz - ? + iho) (2.13) 


X ' 


vho - ch^ Xobo 


For the function oa(z), having a polarity with residue A/Siti, the 
followl 3 ig asymptotic expression is obtained: 


(o.) 


X-*- - <» 


A cos Xq(z - i + 2ihQ) - A cos Xq(z - i) 

iv P~ — ^ (2.14) 

vho - ch^ Xq^O 


Setting A = v(^) d^ and integrating over the contour C 2 _ yields 
the asymptotic expression of the function v(z) = dw/dz: 


dv ^ 

dzJ 


IV 


HC-Xq) exp iXqCz + 2iho) + 


. _» 2 (vho - ch^ Xobo) 

H(Xo) exp [- iXo(z + 2iho)] - H(Xo) exp iXoZ - H(- Xo) exp(- iXoZ^ 

(2.15) 


Finally, from the formula 


, c ^ f dw\ 
y(x) = — Re j ) 

^ g 


it is readily found that for x-^'- sinusoidal waves of length 2jt/Xo 

are formed behind the amplitude of which, after some simple transforma- 
tions, may be represented in the form 


a = 


ch Xoho 


c(vho - ch^ Xobo) 


I H(Xo) exp Xoho - H(- Xo) exp - Xoho | 


(2.16) 
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3. Formulas for Determiniag the Forces 

The forces acting on the contour C are now computed. The lift 
force of the contoixr is denoted hy P, the resistance by and the 
moment of the forces on the contour about the origin by M. These 
forces will be coniputed by the formulas of Chaplygin- Blasius: 


P - iE 



v2q(z) dZ;, M = Re 




zv2q(z) 


dz 


(3.1) 


where C 2 is an arbitrary contour, situated in the region 0>y> -hQ 
and containing the contour C; and Vq(z) is the complex velocity in 

the relative motion obtained by supeiposing on the absolute flow a 
uniform motion of the fluid with velocity c in the direction of the 
negative x-axls. Thus, 


Vo(z) » Vj^(z) + VgCz) - c 

where the contour Cj_ is chosen to lie between C and Cg. 

Formulas (3.l) do not taJse into accoimt the buoyancy 'force of 
Archimedes, equal to gpS, and it moment, equal to -gpSXg, where S 

is the area that bounds the contour C, and x^ is the abscissa of the 
center of gravity of this area. 


The following integral is now conputed; 



vo^(z) dz = / V-l 2 ( 2 ) 



(VgCz) - c)2 dz + 2 / Vj_(V 2 - c) dz 

l/co 


But the first and second Integrals on the right are equal to zero 
because the function Vj_(z) is holomoiphlc outside the contour Cg 

and has at infinity a zero of at least the first order, while the func- 
tion V 2 (z) is holomorphic within the contour Cg. Hence, 


J = 



dz 



(Vi + Vg) dz 




v(z) dz 
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The velocity clreulation about any contour that contains the con- 
tour C is denoted by F , so that 


r = / v(z) dz 
'C 


therefore 


P - iE = - p / V]V2 


(3,2) 


By the use of expressions (2.2) and (2.U), the following expression is 
obtained 



Vi(z) V2(z) dz = ^ 


I r i /[s(- 

vJcg ^Ci L-Uo 


\) exp iX(z + 21hQ) + 


> H(- X) iX(z + 2iho) - H(X) exp[- iX(z + 2iho)] 

— — exp(- Xho) 


V sh XIq - X ch XhQ 


1(- Xq) exp'iXo(z + 2iho) + 1 (Xq) exp [- xXq(z + 2iho)] 
2 vhQ - ch Xoho 

„iv H(Xq) exp IXqz + h(- Xq) exp(- iXpz) 

2 vhQ - ch^ XqIo 



Since the point which belongs to the contour lies within 

the contour with an Interchange in the order of integration and 

by the following formula^ 


1 

2jtl 


r e±±>^ z 

/ 

dc2 


3±iX^ 


There is obtained 
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Vi(z) ¥ 2 ( 2 ) dz » 




|h(- X ) exp(- 2X1 iq) + 


V +X . ^ ^ |H(-X)|^ e^(- 2 x 11 0 ) - i h(X)| 2 ejg) 2X1 iq”| 

“2“ vShXho- \ch\ho J ■ 

nlv H(- Xo)|2 ejj>(- ZX^ho) + I hCXq) P - 2H(X(,) H<-X(,) 


vhQ - cn2 Xobo 

Hence, formula (2.3) assumes the form 


(3.3) 


too 

P - IE = per - X)!^ exp(- 2XhQ) + 

, , , , 1h(-X)| 2 exp(- 2Xho) - 1 h(X)| 2 exp 2Xhol 

(v+X) e^(-Xho) ^ rih - Xh ^ ~ X - ^Xbo l 


XV p 


|H(-Xo)l^ exp(- 2 Xq1io) + IsCXq) exp 2XohQ - 2 H(Xq) H(- Xq) 


vh<3 - ch^ Xoho 


(3.4) 


Separating the real and imaginary parts and adding to P the 
Archimedes force, not taken into account by the Chaplygin-Blasius for- 
mula, results in 


p = per - 


‘i'll 


H(- X)l^,exp(- 2Xbo) + (v + \) exp(- Xho) ^ 


1h(- X)P exp(- 2X1 iq) -|h(\)P exp 2X1 iq 
2(v sh XliQ - X ch XIiq) 


dX + vp 


Im{H(Xo) H(- Xq)} 
2(vhQ - ch^ XQhQ) 


+ gpS 


( 3 . 5 ) 


pv |h(Xq)| 2 exp 2Xoho + |h(-Xq)|^ exp(- aXpUp) - 2 Re{H(Xo) H(- Xq)} 

^ V ho - ch2 Xoho 


( 3 . 6 ) 
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Formula (3.6) may be given another form, namely 


R 


pv ^0^ ®^(“ ^0^0^ I 

^ _ TF7“7 


(3,7) 


It can be readily shown that the total resistance of the underwater 
wing consists only of the wave resistance. In fact, by the following 
well-known formula for computing the wave resistance in the case of a 
fluid of finite depth, 

E = i pga2 (3-8) 

and with the value of the amplitude a from formula (2.16), for- 
mula ( 3 . 7 ) is obtained after some transformations. 

The moment of the acting forces on the contour C is now com- 
puted. When the moment of the Archimedes force is taken into account, 


M = - gpSxc 


+ Re I 



Vi(z) + V2(z) - c 


dz 


(3.9) 


This expression is coniputed in an entirely similar manner to the 
computation of the expression P - iR. 

For very large absolute values of z the following expansion can 
be employed 



Further, 
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z(V 3 - c) dz = 0 


and therefore^ 


M = - gpSx^ + Rep / z Vq_(V 2 “ c) dz 


or^ since the function VgCz) is holomorphic within Ihe contour C 2 


M = - gpSX(, - pc Re i zv(z) dz + p Re / z V 3 _(z) ¥ 2 ( 2 ) dz 




It is noted that 




H*(X) = -^ = - i / exp(- iX^) d^ 


The integrals in formula (3.10) are coit 5 )uted in the same manner as in 
the expression ^3.3)^ and as a resiilt there is obtained the formula 


M = - gpSX(, - pc Re [^lH’(0)]j + p Re 




(- x) H(- x) exp(- SXqIiq) + 


H'(- X) h(x) - H(- X) H'(X)^J ax - 

— fn'(- Xq) h(- Xq) exp(- RXoho) " H‘(Xq) H(Xo) exp SXqRq - 

4(vho - ch2 XqUo) \ 

E'(- Xq) H(Xq) + H(- Xq) H'(Xo))”^ 


E'(- Xq) H(Xq) + B(- Xq) H'(Xq) 
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Formulas (5.5)^ (3.7), and (3.1l) in the limiting case for hQ-*-® 
agree with the formulas obtained by Kochln in reference 2. 

The function H(X) in formulas (3.5), (3.7), and (3,H) does not 
depend on the contour C2_, and for example, the contoTir C or some 

other contour which contains the contour C may be taken for the con- 
to-ux of integration. Moreover, the value of the_function H(X) does 
not change if. Instead of the complex velocity v(z) _of the absolute 
motion, the complex velocity of the relative motion Vq(z) is taken, 

because these two functions differ by a constant c. The properties 
of the function H(X) will be used in the following section. 


4. Examples 

In the preceding sections expressions were found in terms of the 
function H(X) of a numiber of important magnitudes, namely, the anpll- 
tude of the waves formed, the wave resistance, the lift force, and the 
moment of the forces acting on the contour. Thus, the f inaction 


H(X) = / v(^) e::q3(- iX^ ~ j dw exp(- iX^) (4.l) 

l/c Jc 

plays a fundamental part for the problem under consideration. In order 
to compute this function, it is necessary to know the expression for 
the complex velocity, i.e., the solution of the hydrodynamic problem. 

In case the relative depth of the submerged contour C is sufficiently 
large, hovrever, a good approximation is obtained if, in place of the 
function v(z), there is substituted in formula (4.l) the expression 
of the complex velocity which corresponds to the motion of the con- 
tour C in an Infinite fluid. 

Several examples of such an approximate solution of the problem 
will be considered 

1. The motion of a circular cylinder . - The circular cylinder of 
radius b, situated at the depth h under the free surface of the 
fluid, is assumed to move with constant horizontal forward velocity c, 
since the circulation about the contour of the cylinder has a given 
value r. In this case, the characteristic function for the infinite 
fluid is known; 


cb^ + r 
z + hi 27ti 


w (z) = - 


ln(z + ih) 
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Hence ^ 


r(z) =;= 


cb' 


(z + hi)2 2ni(z + ib) 

By formula (4.l) the function H(X) is now constructed: 


(4.2) 


H(X) = 


cb 


(z + hi)2 2ni(z + hi) 


exp - iXz dz 


Since the contour C contains one singulax point z = -ih, there 
is obtained by the theorem on residues 


H(X) = (r + 2itcb^X) exp - Xh 


(4.3) 


With the use of formiila (3.7), the expression for the wave resist- 
ance of the cylinder is obtained 


R = pv 


r sh XQ(hQ - h) + 2jtcXQb^ ch XQ(hQ - h) ^ 


ch^ Xoho - vho 


(4.4) 


and by the use of formula (3.5) the expression for the lift force of 
the cylinder is obtained 


p 


pcF - 


4n(hQ - h) 


pcb^r 
2(ho - h)2 


2(ho - h)3 



+ X) exp(- XIiq) 


(r^ + sh 2X(ho - h) + 4jict)2rX ch 2X(hQ - h) 

vsh Xhg - X ch XhQ 


ax + gps 


( 4 . 5 ) 


The integral component of this formula may be computed by the 
method of mechanical quadratures. In the limiting cases v = 0 
V = this component can be very accurately computed. Moreover, if 
this integral component is considered as a function of the parameter 
a = l/(i/ho) = c^/(gjiQ), it can be shown that for a = 1 this component 

suffers a discontinuity. In the particular case when the radius b of 
the cylinder is taken equal to zero, i.e., when the motion of a vortex 
under a free surface is considered, formulas (4.4) and (4.5) lead to 
the expressions established by Tikhonov. It is noted further that for- 
mulas (4.4) and (4.5) have been derived on the assumption that c2< ghQ. 

For c^>ghQ, no free waves are formed behind the cylinder and the wave 
resistance R is equal to zero. 
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For the moment of the forces exerted hy the fluid on the cylinder^ 
the following expression is obtained hy formula (S.U) : 

PV g*(" ^q) H(" ^q) exp(- aXpUp) - H»(Xq) H(Xq) CXpho ^ 

^ vhQ - ch^ X^hQ 

PV B»(- Xq) H(Xq) - H(- Xq) H«(Xq) 

^ vho - ch^ Xoho 


But from equation (4.3)^ it is evident that 

H* (X) = - hH(X) + 2jtcb^ exp(- Xh) 

H* (- X) = hH(- X) - 2jtch^ ejcp Xh 
Hence, after simple transfoimatlons, 

r sh^ XQ(hQ - h) + jtcb^XQ sh 2XQ(hQ - h) 


M = hR - 2«pcb^v 


ch^ Xoho - vho 


(4.6) 


The point of Intersection with the y-axls of the resultant force 
on the body is determined by the formula 


yo = - 


M 

R 


h + 


2jtcb^ 

r + 27tcb^O cth Xo(ho - li) 


(4.7) 


It is evident that for R>0 this resultant never passes throu^ 
the center of the cylinder. 

2. Motion of an elliptic cylinder . - An ellipse, having a center 
at the depth h and having axes 2a and 2j3 directed parallel to 
the axes of coordinates x and y, is allowed to move with a constant 
velocity c in the direction of the x-axls. The circiilation F is, 
for simplicity, taken equal to zero. In this case, the flow of an 
infinite fluid about the contour C is determined with the aid of an 
auxiliary variable and the formula 
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where r = ^{a, + p)/(a ~ p) and f ul = r is the equation of the cir- 
cle in the u-plane which corresponds to the contour of the ellipse C. 
The exterior of this circle corresponds to the exterior of the ellipse. 
The following function is set up: 



When the substitution u = iv is made, there is obtained 



But by the theory of Bessel functions it is known that 



From the formula 


J_l(z) = - Ji(z) 

and the value of r, the following expression is obtained 

H(X) = 2-acP> exp(- Xh) (4=‘8) 

The confutation is restricted to the wave resistance. By 
formula (3.7), 



18 


MCA TM 1345 


R = 4it^pg3 


2 g + 3 

g “ P ch^ \)liO “ 




(4.9) 


From •fchis foimila, it fpUows that for certain \q and, therefore 
for a certain velocity c <A/gho^ wave resistance is equal to zero; 

i,e., the amplitude of the waves fonned behind the moving body becomes 
zero. This will be the case if the following relation is satisfied: 

Xq f\Ja^ - = sjj. (k=l,2, . . .) 

where Sj^ is the positive root of the Bessel fimction Jj_(s). The 
first root of this function is 


Sj_ = 3.832 


Since the parameter v = g/c^ is connected with Xq by the 
equation 


th XQho 



the first velocity at which the wave resistance becomes zero is deter- 
mined by the formula 


c 


0.51 



p2 th 


3.832ho 



(4.10) 


Moreover, 


3 . 832hf) 

th — — ~ < 1 
A/x^ - 


hence, 


c < 0.51 '\lgNa^ - (4,11) 

In a similar manner a nimiber of other examples may be considered. 
Moreover, as in reference 2, it is possible in this case to set up a 
functional equation for determining the function H(X) and the values 
of the circulation f from the condition of the finite velocity at the 
sharp edge. These equatipns may be obtained by the same method. Their 
final form will be somewhat more complicated as compared with the case 
of the infinite fluid. 
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